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WIMPless dark matter [3| offers an attractive framework in which dark matter can be very light. 
We investigate the implications of such scenarios on invisible decays of bottomonium states for dark 
matter with a mass less than around 5 GeV. We relate these decays to measurements of nucleon- 
dark matter elastic scattering. We also investigate the effect that a coupling to s quarks has on 
flavor changing b s processes involving missing energy. 
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(D . I. INTRODUCTION 

' Numerous cosmological observations require the existence of nonluminous matter that has not yet been observed 
in the laboratory - dark matter. Efforts to explain electroweak symmetry breaking typically introduce as yet unseen 
I— i. weakly-interacting massive particles (WIMPs). Scenarios which economically relate the dark matter and electroweak 
' symmetry breaking problems by positing a WIMP as comprising dark matter are quite common. However, it is 
worth noting that dark matter may not be connected to the weak scale. WIMPless models ^ of dark matter offer a 
natural scenario in which dark matter can have a mass that is not that of typical WIMP. These models incorporate 
O , the minimal supersymmetric standard model (MSSM) with gauge mediated supersymmetry breaking (GMSB) (for 
7^ , ' a review of GMSB, see [2]) and involve a hidden sector that contains the dark matter. This hidden sector dark 
[ matter is unconnected to the weak scale and may have a mass that is quite different in magnitude. In particular, 
PsJ ■ dark matter that has GeV scale mass can be accommodated. Dark matter of this scale has been invoked to explain 
> ' several tantalizing experimental results, notably those of DAMA/LIBRA and SPI/INTEGRAL [lllli]. Sub-GeV 
scale bosons have also been proposed to explain the ATIC Q and PAMELA 0] cosmic ray spectra by providing a 
Sommerfeld enhancement so that the requisite dark matter annihilation cross sections can agree with relic density 
constraints and to suppress decays to hadronic final states |^]. In this note, we examine the implications that a 
particular WIMPless dark matter model has on the decays of bb bound states and on direct detection of nucleon-dark 
' matter elastic scattering. 

WIMPless dark matter is introduced in Sec. II along with a particular model. Invisible decays of bottomonium 
I states in this model are considered in Sec. III. The implications for dark matter-nucleon scattering and its relation to 
• • invisible bottomonium decays are discussed in Sec. IV. In Sec. V, we consider flavor changing b-s transitions in this 
^ model. In Sec. VI we conclude. 

C3 II. THE MODEL 

_ _ I 

We consider a model introduced in [l[ and briefly review it here. We have three sectors: the MSSM, a hidden 
sector, and a so-called secluded sector which breaks supersymmetry (SUSY). SUSY breaking in the secluded sector is 
transmitted to the MSSM and hidden sector with GMSB. In the secluded sector, a chiral field S attains an expectation 
value, {S) = M + 9^F which breaks SUSY. This breaking is transmitted to the MSSM by a messenger field, $, which 
transforms nontrivially under the SM gauge group and interacts with S through the superpotential W = h^S^. The 
spinor components of $ form Dirac fermions with mass Mm and the scalar components have squared masses ± Fm 
with Mm = hM and Fm = hF. This generates MSSM superpartner masses of order 

^ 167r2 Mm 167r2 M ^ ' 
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where g is the largest relevant gauge coupling. In general, a stable thermal relic's density is set by its mass and the 
coupling strength of its annihilation channels, namely, f2 oc rr? j . Combining this and Eq.[l] fixes the density of a 
stable thermal relic in terms of the components of (S"), cx m? / g'^ ~ F^/(167r^M)^. The WIMP "miracle" is that 
we get a relic density of the right order (fi ^ 0.25) if we use weak scale masses and couplings for m and g. 

We now repeat this process with a hidden sector. We introduce a hidden sector messenger, <I>x, coupled to S 
through Wx — hx^xS^x- As before, the spinor components have mass MmX and the scalars have squared masses 
M"^-^ ± FmX with MmX = hxM and F^x = hxF. This sets the scale of hidden sector masses, 

gx Fmx ^ g\ F 
^ 167r2 MmX 167r2 M ^ ' 

with gx the largest relevant gauge coupling in the hidden sector. Then hidden sector thermal relics have a density 
given by Q.x (x m^/gj^ ~ F'^ / (Wtt'^ M)'^ . Thus, we find that a stable hidden thermal relic will have roughly the right 
density, fix ^ 0.25, since its value is set by the ratio F/M as in the WIMP case - a WIMPless "miracle." 

A concrete WIMPless model, with a large hidden sector, was shown in detail in [l^l to be able to reproduce the 
observed dark matter relic density, illustrating that the general formulation above works in practice. 

In the standard GMSB scenario, there is a problem with gravitino production in the decay of the dark matter 
candidate due to the lightness of stable Standard Model (SM) particles. This lightness is due to extremely suppressed 
Yukawa couplings and, as in we assume that such a situation does not exist in the hidden sector. 

It is phenomenologically interesting to consider a case in which the SM is coupled to the hidden sector dark matter 
through some connector. These interactions often occur in intersecting brane models in which a connector has both 
hidden sector and SM gauge quantum numbers [ll|. We consider a case in which the hidden sector dark matter is a 
scalar, denoted by X. It is coupled to SM fermions, written as / here, through a chiral fermion, Y , via the interaction 
Lagrangian, 

Cint=XfXYLfL + XfXYRfR . (3) 

This model allows for dark matter masses in the range 10 MeV < mx < 10 TeV. The lower bound is set by the 
requirement that X be in thermal equilibrium during freeze out so that the expression for the relic density of X, 
fix cc m\, is valid. The upper bound is found by enforcing that gx remain perturbative while giving roughly the right 
relic density. We are interested here in very light dark matter, mjc < 5 GeV. The connector Y is subject to limits on 
extra colored particles from the Tevatron: my > 258 GeV [I4I . Perturbativity requires my < 500 GeV since Y gets 
contributions to its mass from electroweak symmetry breaking. The relic density estimate for X, fix oc F^/ (IQn'^M)'^, 
is not affected as long as A/ < (?wcak- This Lagrangian generates several interesting signals which we describe below. 

III. INVISIBLE xbo DECAYS 

The Lagrangian of Eq. [3] induces the reaction bb XX when we take f = b. Similar interactions were considered 
in [1^ involving J^*^ = ^ , 1 quarkonium states. In this case the relevant quarkonium state has J^'^ = 0"*"^. In 
the limit of large my the t-channcl (shown in Fig. [T]) contribution to the amplitude is 

iMtibb ^ XX) ~ ^^v{p')u{p) . (4) 
my 

In this limit, the M-channel contribution is the same and so the total amplitude is given by 

iM{bb -> XX) ~ -^v{p')uip) . (5) 
my 

In dealing with nonrelativistic bound states it is useful to write the quark spinors in a chiral basis to linear order in 
their three-momenta, 
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FIG. 1: f-channel bb XX diagram that leads to xto XX. 



where p is the center-of-mass momentum of the b quark, E is its energy, and f and ^' are two component spinors of 
the b and b quarks, respectively. Then, 



iMibb^ XX) 



where 



my 



X = Tr 



my 



P ■ X 



V2 



(8) 



(9) 



describes a spin-1 configuration of the quarks. Thus, this matrix element induces a coupling of XX to the '^Pq bb 
states, the X6o(l-P) and Xfco(2P). These states' wave functions are given by 



V3 



r \ 4tt 



Rnp{r) 



Properly accounting for the bound state wave function, the matrix element for the Xbo to decay to XX is 



MixbainP)^ XX) = 



XbO 

8X1 



8X1 



y/3M^my 
12i 



-3z^/— Kp(0) 



Thus the decay rate is 



TixwinP) ^ XX) = 



2Af. 



Xba 
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\Mixbo^XX)\' 



27 \Kpm 

27r2 Ml^m^ 



with = 1 — Am'j^ l-^Xbo ■ Using the expressions for the widths of the Xbo states 'l4| , 
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we can find the branching ratios to a pair of X bosons, 

B{xbQ{nP) ^ XX) ~ 1.3 X 10-4^2 ^ 

with 

\ my J 

XboS are produced in radiative decays of the T{2S) and T(3S'), T{2S) — > 7Xbo(l-P) and T(3S') 
Ts win decay to ^XX through the Xbo and we find the branching ratios 

B{T{2S) ^ jXX) = B{T{2S) ^ jXbo{lPmxbo{lP) -> XX) 
~ (4.9±0.5) X 10"^^^ , 

B{T{3S) ^ -/XX) = B{Ti3S) ^ 7Xw(2F))6(xbo(2F) ^ XX) 
~ (7.7±0.8) X 10"^^^ , 



(20) 



(21) 



7X60 (2P). Tlius the 

(22) 
(23) 
(24) 
(25) 



where the error shown is only due to that of the measured branching ratios, B{T{2S) IXboi^P)) = (3.8 ± 0.4) % 
and S(T(35) 7Xbo(2P)) = (5.9 ± 0.6) % ^ 



A major background for this search at e'^e colliders is radiative Bhabha scattering, e^e 
electron and positron go undetected. The cross section is estimated to be (Eq. IA.36P 



7e^e 



da 



dcos 9 



log 



log 



S2 (1 - c) 



2,771^ ^min ^max 



(1 



where the 



(26) 



where ^/s is the center-of-mass energy and m is the mass of the electron, c is determined by the requirement that 
the electron and positron are unseen: |cos0+__| > c where 9+ (9-) is the angle the positron (electron) forms with the 
beam. The photon has energy between fcmin and fcmax and 9 is the photon's angle with the beam direction. 

The photon in the T(2S') JXboi^P) transition will have an energy of 160 MeV while in the T(3S') jXboC^P) 
case it will be 123 MeV. We have in mind the CLEO experiment and assume a resolution on the photon's energy of 
3.5%, and that c = 0.9 for the electron and positron to be unseen. The cross section at both the T(2S') and T(3S') 
resonances is then 



da 

dcos I 



(75 pb) 



1 + cos^ 



(1 - C0S2 6')2 



(27) 



The CLEO experiment has collected 1.3 fb^^ on the T{2S) resonance and 1.4 fb ^ on the T(3S'), producing 9.32 x 10^ 



T(2S')s and 5.88 x 10^ T(3S')s 46]. The signal for T -> -fXbo ^ l^X wiU be distributed as (3/8)(l + cos^ 



We 



require |cos^| < 1/V2 to cut down on the radiative Bhabha background which is extremely peaked along the beam 
direction. We expect this sample to contain 



T(2S') JXboilP) iXX events and 



(29 ± 3) 



(28 ±3) ^2 



(28) 



(29) 



T(3S') jXboi'^P) jXX events. Setting these equal to the variation in the background due to radiative Bhabha 
scattering from Eq. [27| implies that one can probe 



C > 4.3 

on the T(2S') resonance and 

C>4.4 

on the T(3S') resonance which correspond to the branching ratios 

B{T{2S) -yXX) > (9.1 ± 0.5) x 10~^ 
B{T{3S) jXX) > (15 ± 0.5) x 10"^ 



(30) 
(31) 

(32) 
(33) 
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FIG. 2: Diagram that leads to spin independent X-proton cross section. 



If TOy = 400 GeV and mx = 1 GeV these limits on ^ imply a limit can be set on the coupling of Y and X to b quarks 
of Ab > 2.1. This is a little larger than the maximum value of Xb that would not disrupt the relic density estimate of 
X. We have assumed here an efficiency of 100%. The limits on the branching ratios should scale roughly with the 
inverse of the square root of the efficiency. 

Regarding the reach of other e+e^ collider experiments in setting limits on B{T •^XX), we note that the BaBar 
experiment collected over 100 X 10^ T(35) decays [13 • The reach in the branching ratio, given similar resolutions 
and efficiencies, scales as iV~^/^, with N the number of decays, so one would expect BaBar to be able to probe 
B{T{3S) -fXX) - 10"^. A Super-B factory like that envisioned by the BELLE collaboration with 10 - 20 ab"^ 
collected on the T(4S') could collect 10^ T(3S')s and limit B{T{3S) jXX) - 5 x 10"^ 

IV. NUCLEON-X SCATTERING 

The Lagrangian of Eq. [3] can also lead to nucleon-X interactions like the one shown in Fig. O This leads to a spin 
independent cross section [l^ 



m%[ZBl + iA-Z)B-f 
47r ^2 (^jjij^ 4_ mx)^ {rriY — mx)'^ 



■ -b n ^ ' \ / OJ ('Q/|^ 



where A is the atomic mass of the nucleon, Z is its atomic number, and B^'" = {p,n\bb\p,n). This matrix element 
has been calculated using the conformal anomaly [l^, and we write _B^'" = (2/27)m7v/^'"/mb and, as in [l^, we take 
/I = /g ~ 0.8 = fg. Then the cross section becomes 

An \27 J ml{mN + mx)^{iTiY — mxY 



-si-^ ^..^ . (35) 



Setting a limit on B{T{2S) JXbo ~^ iXX), through Af,, corresponds to setting a limit on crsi. We have plotted in 
Fig.[3]the reach in ctsi that can be obtained for particular limits on B{T{2S) JXbo jXX). We use a messenger of 
mass my = 400 GeV. The CRESST experiment [l^l sets the most stringent limits on csi for rnx < 5 GeV. We have 
also indicated the range of asi and mx that are compatible with the dark matter interpretation of the DAMA/LIBRA 
data if one allows for dark matter streams Q. However, the reach in ctsi for a given limit on S(T(2S') 7X6o —> jXX) 
is nearly independent of the value of my. This is because both crgi and ;B(T(25) JXbo iXX), for large my, 
depend only on the quantity A^/rriy. 

Fig. [3] and the limits from Eqs. [32land [33] show that current CLEO data, in setting a limit on the branching ratio 
B{T — > jXX) > lO"**, can complement data from direct detection experiments for dark matter particles with a mass 
mx < 1 GeV. A limit B{T jXX) > 10^^ could be possible at a Super-B factory and would be quite competitive 
with those from direct detection. 

V. b-s TRANSITIONS 

The Lagrangian of Eq. [3]can give rise to an effective b-s quark interaction if both A^ and As are nonzero. If so, to 
lowest order in my, one finds an effective interaction term, 

C,f, = ^^sbX^ . (36) 
my 
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FIG. 3: Spin independent X-proton cross section as a function of mx with contours of constant B(T(25) 7X60 — > yXX) 
indicated. The CRESST limit [20l ] is shown in red (solid) while the area between the red (solid) and ereen (dashed) lines is 
consistent with the dark matter interpretation of DAMA/LIBRA if dark matter streams are included [3]. 
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FIG. 4: Diagram relevant for B+ -> K+XX. 



This will induce the decay K~^XX through the diagram shown in Fig. [4l Using this process to constrain other 

models of dark matter has been considered in '21'| . The amplitude for this decay, mediated by the effective Lagrangian 
of Eq. Ea is 



M(B+^K+XX)^^^^^^^^fo{q') , 
my rrib — nis 

where = {pb — Pk)"^ ■ /o(<Z^) is estimated, using light-cone sum rules, to be [l^ 



/o(g') = 0.3 Exp 



0-63 ( 4tt] - 0-095 ( 4w 



0.591 



.2 \ 3 



Ml 



The differential decay rate is then 
dV 



, iB+ K+XX) = — — ^ 



niY 



' ( Ml~Ml 
\ rrib-ms 



fl{q^)X'/\q\Ml^Ml)Jl 



4m 



X 



with 



Ai/2(q2^ ^^2 ^ ^ ^ ^,^4 _^ ^,^4 _ 2q^Ml - 2q^Ml - 2Af|M 



1/2 
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FIG. 5: The upper limit on \p\, defined in Eq.|42l as a function of mx from B {B+ K+vv) < 1.4 x 10"^ [H] and from the 
contribution to Arris from X and Y exchange being less that the measured value of 17.77 ± 0.12 ps~^ for my = 400 GeV. 




FIG. 6: One diagram that gives a contribution to Arris. 



Using mi, — 4.2 GeV and iris — 100 MeV, we integrate over and find the branching ratio, 

B{B+ ^ K+ XX) = {1.0x10'^) \p\^ F{mx) , 

where we define 

/ 400 GeV - 

and F{mx) describes the allowed phase space as a function of mx, 



(41) 



(42) 



F(mx) 



(Mb-MkT 



J^''--''-^ dq^P,{q^)\y^{q\Ml,Ml) 



(43) 



Using the Umit B {B^ K^vv) < 1.4 x 10~^ [1^, we can set an upper limit on \p\. For mx < 2 GeV, we find 
p < 10^''. The upper Umit on p as a function of mx is shown in Fig. O 

There is also a contribution to Eg — Eg mixing from the Lagrangian of Eq. [3] at one-loop through bs — > bs. One 
diagram that contributes is shown in Fig. [SI The relevant portion of the amplitude for this is 



M (bs bs) 



Z2T:'^m\ 



[5(y) + /i(y)]r7W7^s" + 



9{V) 



7m7 



(44) 
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where a, P are color indices, y = m^/my and 



9iy) = 

h{y) = 



(1 - yf 

2(y-l)-(l + y)log(y) 



(i-yr 

In the limit that mx <C my, using the vacuum saturation approximation, this yields 



{B,\M\Bs) 



m 



mi, + rris 



which gives a contribution to the mass difference 



Am, 



log 



2887r2m2, 



X 



8 + 5 



/ Ah 



(1.82 ± 0.12 X 10^ ps^i) n - 0.08 log 



mx \2 /400 GeVV 
1 GeV/ V j 



(45) 
(46) 

(47) 

(48) 
(49) 



where we have used /b^ ~ 231±15 MeV [2J|. Conservatively using the total measured value of Am^ = 17.77±0.12 ps^^ 
[25| as an upper limit on the expression above, we can set an upper limit on ignoring any phase difference between 
Ab and Ag. This limit, for my = 400 GeV, is shown in Fig.O For values of mx where — > K^XX is kinematically 
allowed, its constraint on p is much stronger than that of Awg. 

Limits from h-s transitions give stricter constraints on couplings of the form of those in Eq. [3] than those involving 
only couplings to h quarks. We are led to speculate that there may be a Cabibbo-like suppression of the coupling of 
a heavy chiral fermion, Y ^ to scalar dark matter, X, and an s quark relative to that involving a 6 quark. 



VI. CONCLUSIONS 



Bosons with a mass on the GeV scale that are weakly coupled to the SM have been implicated in a number of 
experimental observations, and may even comprise the dark matter. WIMPless models of dark matter can naturally 
accommodate scalar dark matter with a mass in the GeV range. Scenarios with light dark matter can have many 
implications for decays of mesons that involve undetected final states (missing energy). We have seen here an example 
of a simple WIMPless model involving a fermionic messenger and scalar dark matter that gives novel signatures. The 
conclusions drawn here, however, apply to any model involving an interaction of the form of that in Eq. [31 Existing 
constraints from collider experiments can verify those from direct detection experiments or, in the case of very light 
dark matter, be complementary to them. Comparison of effects involving flavor changing b-s transitions and those 
only involving decay of hh states indicates that it is likely that dark matter in this scenario is preferentially coupled 
to h quarks. In this case, studies of the decays of bottomonium states involving invisible final state particles provide 
the opportunity to test these models. 
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APPENDIX: CALCULATION OF THE RADIATIVE BHABHA CROSS SECTION 

The squared amplitude for radiative Bhabha scattering, in the limit that the electron mass m = 0, is (Eq. (2.60) 
of [2^) 

\M\^ = - "^^ ~ ""^^^ [ss' {s^ + s'^) + W + t'^) + uu' {u" + u'^)] , (A.l) 
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where 



s = {p++P-f , s' = {q+ + q-f , 
t = ip+- q+f , t' = {p- - q-f , 
U = {P+- Q-f , u' = {p- - q+f , 



1+ 



q+ ■ k q- ■ k 



Vr, = 



P+ 



p+ ■ k p- ■ k 



(A.2) 
(A.3) 
(A.4) 
(A.5) 



We take p_ as the four- momenta of the incoming electron (positron), q_ as the four-momenta of the outgoing 
electron (positron) , and k as the photon's four-momentum. We write the four-momenta in the center-of-mass frame, 
neglecting the electron's mass, as 



P+ = 



g+=g+(l,n+) , q-=q°_{l,n-) , 

h± = sin 6± cos (l)± x + sin 9± sin (l)± y + cos 9± z , 



k 



(l,n) , h = sm6cos(f) X + smOsiiKp y + cos6 z 



We find 



l-x 



1- (a;/2)(l - n-n+) 
l-{x/2){2-x){l-h-n+) 



1- {x/2){l-h-n+) 

(1 — x)n+ + x{l — {x/2){l — n ■ n+))n 



l-{x/2){2-x){l-h-h+) 



(A.6) 

(A.7) 
(A.8) 

(A.9) 



(A.10) 
(A.11) 
(A.12) 



We are interested in cases where the photon's energy is small compared to the center-of-mass energy and the electron 
and positron are slightly deflected. That is, we are interested in cases where x is small. To first order in x, 



l--(l+n-n+) 

l-|(l-n-n+) 



{1 — xh ■ n+)n-)- — xfi 



We use these to find Vp and Vg to lowest order in .x. 



cos 9 



xws / sm- 



x^/s J sin^ 6 



4 

X^/s 



1 — (n • n+)^ 



1 - (n • n+f 



(A.13) 

(A.14) 
(A.15) 

(A.16) 
(A.17) 
(A.18) 
(A.19) 



The dominant contribution to the cross section comes when cos^+ ~ 1. We write cos^+ = 1 — and expand Vg 
to lowest order in a, 



x^s 
4 



cos ^ -I- a sin ^ cos (f)+{l -|- 2cot^ 9) 



sm^9 

z + a(cos 4>+x-\- sin (/>+y + 2 cot 9 cos (/>+2) 
sin^6l 



(A.20) 
(A.21) 
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where, without loss of generaHty, we have set cos0 = 1. Then we find 

'16a^\ rsin2 6'cos2 0+(l + 2cot2 6')2 - 1 - 4 cot^ 6* cos^ 



sm 



16a" 



1 — sin 9 cos^ 



sm 



(A.22) 
(A.23) 



Now, looking at the other Lorentz invariants that make up the amplitude we also find to lowest order in x and a, 

s' ~ s (A.24) 



(1 — cos ( 



sa 



Collecting all of this we find the squared amplitude to lowest order in x and a. 



\M\' 



The differential cross section is given by 



da ~ 



2i0g6\ r 1 _ sin^ 6) cos^ 



sm 



(27r)5 8s3/2 



(A.25) 
(A.26) 

(A.27) 
(A.28) 



where (p^, is the momentum of the final positron in the final electron-positron rest frame, (k, 17^) is the momen- 
tum of the photon in the center-of-mass frame, and ~ s' — s — 2y/sko. To lowest order in x, the center-of-mass 

frame and final electron-positron frame do not differ so we will carry out the phase space integral in the center-of- 
mass (i.e. we expand the Jacobian to zeroth order in x which is simply unity). Then, to lowest order in x and a, the 
differential cross section is 



da 



1 



2i°e'^ 
(27r)5 8s3/2 \sx^a^ 



1 — sin^ 6 cos^ ( 



Vs \ / ^/sx \ / y/j dx 



dQ^dflk 



(27r)5 



24e6 



1 — sin 6 cos^ (f)+ 



dx 



X\/l — X 



(A.29) 
(A.30) 



Now we write dft+ — dcos6'_|_d0+ = {2/ s)dtd(t)+. We also express in terms of t. Integrating over azimuthal angles 
we find 



da 
dcos9 



1 



(27r)3 



1 - i sin^ ( 
sin'' 6* 



1 + cos^ 9 
(l-cos2 6l)2 



dx dt 



X\/l ~ X t 

dx dt 



x^/\ — X t 



(A.31) 
(A.32) 



This will diverge if we use Eq. (IA.25P to find the upper Hmit on the integral over t by setting cos 9+ — i.e. tmax — 0. 
To estimate the limits on t we expand to the next nonvanishing order in x and include the electron mass. We obtain 



Setting (cos0+)j. 



t~-2p'J_g° (1-COS6I 
c and {cos9j^)^^^ = 1 we find 

imin =i - (0 (1 - c) 



{l + fi-n^) cos( 



(1 -f- cosf 



Integrating over t now gives us 



da 



dcos9 



32a^ 



1 + cos^ 



(1 - cos2 9)^ 



log 



2s (1 - c) 



(1 + cos2 9) 



dx 



(A.33) 



(A.34) 



(A.35) 
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We integrate over photon energies now and find 



da 



dcos9 



32a-' 



1 + cos2 e 



(1 - cos2 



log 



log 



(1-c) 



2m2 fcminfcmax (1 + COS^ 6*) _ 



(A.36) 
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